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Boundary  Layer  Ducting  Effects  on  the  Transmission 
of  Infrared  Wavelengths 

R.  Felnberg,  H.  V.  Hltney  and  H.  G.  Hughes 

Electromagnetic  Propagation  Division 
Environmental  Sciences  Department 
Naval  Ocean  Systems  Center 

An  experiment  was  conducted  to  determine  whether,  ducting  phenomena 
(gradients  In  refractlvlty)  occur  close  to  the  ocean  surface  at  infrared 
wavelengths.  A  CO^  laser  beam  (10.6  ym)  was  transmitted  and  received  over 
an  8.5  km  path  In  San  Diego  Bay.  For  the  cases  Investigated  (sea  vrarmer  than 
the  air),  the  results  of  the  study  Indicate  that  for  air  temperature 
gradients  less  than  l°C/meter,  the  distance  to  the  Infrared  horizon  was  only 
60%  of  the  normal  geometrical  horizon. 
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I.  INTRODOCTION 

Predicting  the  effects  of  changing  weather  conditions  and  atmospheric  constit¬ 
uents  on  sensor  performance^  remains  an  outstanding  Kavy  problem  In  the  design  and 
development  of  improved  electro-optical  (BO)  systems.  Of  special  interest  here 
are  the  foxvsrd-loo)clng  infrared  (flir)  imaging  systems  v4xose  performance  is 
severely  limited  in  regions  of  high  humidity.  This  is  the  case  to  be  expected 
for  fllr  systems  operated  from  periscopes  near  the  ocean  surface. 

SC  systems  performance  close  to  the  ocean  surface  can  also  be  affected  by 

differ€.rce8  between  the  sea  (T  )  and  air  (T  )  temperatures.  Moderate  temperature 

8  a 

gradients  can  result  in  large  gradients  in  the  optical  and  infrared  refractlvlty. 
This  can  give  rise  to  refraction  (or  bending)  of  imaged  rays  whl<di  may  either 
increase  or  reduce  (depending  vqpon  the  sign  of  the  temperature  difference)  the 
detection  ranges  of  fllr  systemc. 

In  the  standerd  atmosphere  (%Aiere  the  temperature  lapse  rate  is  of  the  order 
of  0.6^C/100  metres)  light  rays  are  bent  downward  but  with  a  radius  of  curvature 
greeter  than  that  of  the  earth.  TO  an  observer  on  the  earthy  the  rays  appear  to 
be  bent  slightly  upward.  The  standard  optical  horison  Is  therefore  only  about 
8%  greater  than  the  geoerntrlcel  horison.  When  a  target  is  located  over  the 
standard  optical  horison  as  depicted  in  Figure  I-l  (distances  not  drawn  to  scale) , 
it  would  be  outside  the  field  of  view  of  a  flir  system  operating  close  to  the 
ocean  surfsca  and  escape  detection.  When  the  tenperature  gradient  (T  <T  )  is  of 

8  A 

the  order  of  l^C/mstra,  strong  bending  of  the  rays  towards  the  earth  can  occur. 

This  is  Illustrated  in  Figure  T-2  where  infrared  radiation  emitted  by  a  surface 
target  la  shown  ducted  to  a  detection  system  beyond  the  standard  optical  horison. 

When  T  >T  ,  rays  are  strongly  bent  upward,  thereby  shortening  the  optical  horizon 

8  ft 

and  limiting  the  detection  range  as  shown  in  Figure  1-3. 


Figure  i-l.  Target  beyond  etenderd 


optical  horison. 
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Figure  1-2.  Target  within  extended  horizon  (T 
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There  is  adequate  evidence  of  the  occurrence  of  ducting  phenomena  (mirages) 
for  vlsihle  wavelengths.  Examples  are  shown  in  Figure  Z-4  of  how  objects  can  be 
distorted  when  viewed  from  different  heights  above  the  sea  surface.  The  pictures 
(a) ,  (b)  axid  (c)  show  objects  In  San  Diego  Bay  ^en  viewed  from  heights  above 
the  water  surface  of  approximately  1.5,  .3  and  .15  metres,  respectively.  Seen  In 
(a)  and  (b)  are  an  abandoned  floating  house,  a  "saucer”  shaped  house  boat,  and  spar 
and  ball  buoys,  nxe  dletanee  to  the  stilt  house  Is  approximately  4  km  and  the 
saucer  boat  Is  3  km.  At  a  height  of  1.5  metres,  inferior  Images  of  the  stilt 
house  and  saucer  boat  are  beginning  to  form  and  become  more  pronounced  when  the 
viewing  height  is  lowered  to  .3  metre.  Mote  also  the  position  of  the  spar  buoy 
and  Its  image  relative  to  the  stilt  house  in  (a)  and  (b) ,  and  also  the  Inverted 
image  of  the  ball  buoy.  When  vlsfwed  from  a  height  of  .15  metres.  In  (c) ,  an 
apparition  of  the  saucer  boat  loose  above  the  horlson  as  a  two-image  inferior 
mirage  indicating  the  presence  of  a  strong  positive  gradient  in  optical  re  free  t- 
Ivlty  (ducting)  close  to  the  water  surface,  if  similar  ducting  conditions  exist 
fbr  XR  wavelengths,  the  detection  ranges  of  flir  systems  may  be  enhanced.  However, 
the  resulting  mirage  effects  may  also  degrade  the  system's  target  identification 
capabilities.  It  also  must  be  determined  tdiether  the  effects  of  water  vapor 
absorption,  wind,  and  sea  surface  roughness  counteract  the  ducting  conditions. 

Baported  In  this  technical  note  are  the  reaults  of  an  investigation  to 
determine  If  dwtlng  influences  the  Lang  range  propagation  of  ID  wavelengths 
close  to  the  ocean  surface.  Presented  in  the  following  sections  are  the 
analytical  background  of  refraction  and  computer  programs  vdxich  were  developed 
tc  model  the  refractive  effects;  the  experlmentel  apptoeeb;  the  instrumentation 
description;  meesurements  and  data  analyses;  and  conclusions.  The  work  was 
ecconplished  under  MOSC  IR/IEO  funding,  problem  number  ZE625and  covers  the 
period  from  1  October  1977  to  30  Septesber  1978. 


Figure  1-4.  Various  objects  viewed  near  the  sea  surface.  Observer  height  is 
1.5  m  in  (a),  0.3  o  in  (b)  /  emd  0.15  m  in  (c) . 


XX .  BACKGROUND 

A.  AtiiiOBph>rlc  Attenuation 

Th*  atanosphere  attiinuataa  light  tfavaa  by  two  tnechanlems#  scattering  and 
absorption.  At  visible  wavelengths,  extinction  is  primarily  due  to  scattering 
by  atnospliric  gases  and  aerosols.  Infrared  wavelengths  however,  lose  energy 
both  by  aerosol  scattering  and  absorption,  and  by  gaseous  absorption.  The  magni¬ 
tude  of  the  contxxbutlon  to  total  extinction  by  each  mechanism  depends  upon  the 
optical  properties,  atstospheric  concentration  and  temperature  of  the  molecular 
or  particulate  species. 

In  the  8-12  ym  region,  the  dominant  contribution  to  extinction  by  noleoular 
absorption  is  that  from  water  vapor,  eot  high  humidities,  refrectlve  effects 
which  ars  normally  observahle  in  the  vieible.  may  not  be  epparent  when  viewed 
at  IR  wavelengths  due  to  the  high  atten\uition.  The  effects  of  temperature  and 
relative  humidity  on  the  transmittance  (ratio  of  power  received  to  that  trans¬ 
mitted)  of  a  10.6  ym  wavelength  soi'xoe  over  an  8.5  km  path  is  shown  in 
Figure  Zl-1.  The  calculatione  were  performed  using  ths  lOWTRAM  3B  (Selby,  et  al.) 
atmoapherlc  trensmisslon  model  and  aerosol  effects  were  not  considered.  The 
limitations  of  humidity  on  transmittance  ere  readily  apparent.  For  example,  at 
an  air  temperature  of  20°c  end  85%  relative  humidity,  only  10%  of  the  transmitted 
power  reaches  the  detector.  Oonversely.  for  haze  conditions  and  low  relative 
humidities,  the  scattering  by  aerosols  may  mask  ducting  at  visible  wavelengths 
but  leave  tha  infrared  wavelengths  unaffected. 

D.  Refrectivlty 

Light  waves  travel  through  the  atmosphere  with  a  velocity  v  «diich  is 
Isss  than  that  for  free  space,  c,  by  a  factor  of  1/n  inhere  n  is  the  index  of 
refraction.  If  n  changes  along  the  wavefront  (normal  to  the  direction  of  propa- 
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RELATIVE  HUMIDITY.  % 


Figure  II-l.  Transmittance  of  8.5  km  path  at  10,6  ym  vs.  relative  humidity  for 
various  teiRi>erature8 . 
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gatlon) ,  different  portions  of  the  vave  travel  vLth  different  velocities  causing 
Isending  or  distortion  of  the  wavefront.  In  air  however,  light  waves  are  bent 
through  such  small  angles,  that  objects  appear  'listorted  only  when  viewed  from 
distances  of  several  kilometres. 

Becav  -;e  the  index  of  refraction  of  air  is  close  to  unity  (1.00032  at  STP) , 
it  is  convenient  to  define  the  guantlty  N,  given  by 

N  -  (n-1)  •  10®  (1) 

which  ic  the  refractive  modulus  or  refractivity.  For  light  waves,  the  refractivity 
of  air  depends  upon  ^he  air  density  and  not  the  water  vapor  pressure  as  it  does 
for  radio  frequencies.  For  a  wavelength  X,  the  refractivity  is  given  by 

N  =  (n-1)  •  10®  -  77.6  J  +  0.584  (2) 

^  tx^ 

where  F  is  the  barometric  pressure  in  millibars  and  T  is  the  air  temperatvure 
in  degrees  Kelvin.  Wiile  there  is  a  slight  wavelength  dependency  on  N  (3%  dif¬ 
ference  between  N  values  at  .Sum  and  10  u  n) ,  it  is  the  gradients  in  T  and  P  that 
determine  ducting,  and  all  optical  frequencies  should  be  equally  affected.  The 
vertical  gradients  in  T  and  P  can  be  introduced  into  equation  (2)  by  logarithmic 
differentiation  with  respect  to  the  vertical  ordinate,  Zt 

.  i  (3) 

n-1  a  z  p  az  T  az 

Replacing  the  index  of  refraction  by  c/v  and  introducing  the  equation  of  state 
and  the  hydrostatic  equation  into  the  pressure  gradient  term  of  (3)  we  arrive 
at 

C  ^  .  la  +  i  3^  (4) 

v(c-v)  az  RT  T  az 

where  K  is  the  molecular  v/eight  of  air,  g  is  the  acceleration  cf  gravity  and  R 


is  the  universal  gas  constant. 


C.  Apparent  Elevation  for  Linear  Gradients 

The  difference  between  the  true  and  apparent  elevation  of  an  object  can 
be  now  estimated  if  we  assume  the  temperature  gradient  to  vary  linearly  with 
height  (Fleaglcf  1963).  In  Figure  11-2,  assume  the  temperature  to  decrease  with 
altitude  and  the  object  and  observer  to  be  at  the  same  height  above  a  flat  eea±h. 
Since  the  density  decreases  with  altitude «  the  vqpper  portion  of  a  wave  launched 
at  an  angle  from  the  object  will  travel  with  a  higher  velocity  than  the  lower 
portions.  The  ray  (wave  normal)  will  be  bent  downward  with  a  radius  of  curvature, 
r,  such  that  the  c^ject  as  seen  from  a  distance  S  will  appear  to  be  at  an  elevated 
position  h.  If  horizontal  gradients  in  temperature  are  neglected,  and  for  a 
fixed  value  of  S,  the  time  for  the  wave  front  to  travel  from  the  object  to  the 
observer  is  approximately  independent  of  r,  or  dr/v  is  a  constant,  and 


1  da  1 

?  dr  r 


i  52. 

V  dr 


(5) 


For  a  nearly  horizontal  ray,  h  can  be  considered  a  prolongation  of  r  in  the 
vertical  z  direction  such  that  dv/dr  =  dv/dz  and 

C..6  -  =1  - 

2  r 

where  cos 9  has  been  expanded  In  a  Taylor  series  assuming  6<<1  when  h  <<r. 
Therefore, 


(6) 


1. 

r 


2h 

_2 


(7) 


Combining  equations  (5)  and  (7)  and  since  9  varies  only  slightly  with  r  we  have 


(upon  substitution  into  equation  (4)) 


.  ,n-l.  S  ,Ito  .  3T. 


(8) 
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Equation  (8)  e:q>resses  the  height  difference  with  respect  to  a  Cartesian  coordin¬ 
ate  system.  To  relate  h  to  sea  level  the  correction  factor  =  a  (sec  (s/a) -1)  must 
be  subtracted  where  a  is  the  radius  of  the  earth.  This  correction  is  slight  for 
short  distances  but  amounts  to  7.8  cm  at  1  ]aa. 

The  apparent  elevation  of  an  object  vs.  temperature  gradient  for  different 
distances  is  shown  in  Figure  II-3.  The  change  in  elevation  is  less  sensitive 
to  the  actual  temperature  than  to  its  gradient.  This  is  evident  at  a  distance  of 
1  Ion  where  the  values  of  h  for  different  temperatures  are  compeured.  The  figure 
demonstrates  the  effect  of  the  sign  of  the  temperature  gradient;  in  general  positive 
gradients  (temperature  Increase  with  height)  elevate  the  apparent  position  of  the 
dbject  while  for  negative  gradients  it  is  depressed.  The  effect  of  the  pressure 
gradient  and  6  account  for  the  fact  that  the  lines  do  not  pass  through  the  origin. 

D.  Arbitrary  Gradients 

So  far,  we  have  only  concerned  ourselves  with  the  effects  of  a  linearly 
varying  refractivlty  on  the  propagation  of  optical  rays.  However,  in  order  to 
maintain  conplete  generality,  it  is  also  necessary  to  consider  the  refractive 
effects  resulting  from  an  arbitrary  refractivlty  profile. 

In  general,  an  «u:bitrary  refractivlty  profile  may  be  decomposed  into  any 
nTiniber  of  piecewise  linear  segments.  Snell's  law  for  spherical  boundary  sur¬ 
faces  (surfaces  of  constant  index  of  refraction  around  the  earth  are  spheres) 
may  then  be  applied  and  the  ray  bending  at  each  interface  calculated.  So  starting 
with  a  ray  of  loiown  position  and  inclination,  it  is  possible  to  successfully 
raytrace  its  path  through  the  atmosphere.  A  computer  program  called  HAYTRACE 
(Appendix  A)  was  developed  to  perform  these  calculations. 

Figure  II-4  shows  a  double  image  of  a  heliuii\-neon  laser  beam  propagating 
close  to  the  oceam  surface  and  viewed  from  a  distance  of  4.6  Nm.  Using  the 
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Figure  II-3.  Apparent  elevation  of  an  object  vs.  temperature  gradient  for 
various  distances. 
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assvnned  exponential  temperature  profile  shown  In  Figure  IZ-S,  the  oosputer 
program  was  able  to  generate  the  raytrace  shown  in  Figure  I 1-6.  Notice  that 
there  are  two  rays  which  leave  the  source  and  intersect  at  the  observing  dls- 
temce  of  4.6  Ntn.  An  observer  at  this  distance  located  at  19  ft.  above  the 
surface,  vould  project  each  ray  back  along  the  direction  It  reached  him  end 
form  a  separate  image. 

Exeunples  of  Infrared  laser  beam  Images  are  shown  In  Figures  II-7.  The  linages 
were  taken  with  a  cynard,  Inc.  Model  210  Fast  Scan  Infrared  Camera  which  was  borrowed 
from  NOSC  Code  7313.  The  system's  detector  is  sensitive  to  the  8-14 ym  wavelength 
band  and  uses  a  trl-metal  <Hg  Cd  Te)  detector  cooled  by  liquid  nitrogen. 

Figure  II-7(a)  shows  a  double  Image  of  a  00^  laser  located  8.5  km  away.  The 
image  was  obtained  vhen  the  tide  was  coming  in  and  the  detector  was  situated  5.25 
metres  above  the  «rater.  The  water  was  w'armer  than  the  id.r  and  at  this  point,  the 
image  shape  was  changing  rapidly  with  time.  In  Figure  II-7(b),  which  was  taken 
a  fe'..  minutes  later,  the  double  beam  image  had  disappeared.  These  exanples  Illus¬ 
trate  that  image  distortion  can  occur  at  IR  wavelengths  as  well  as  in  visible 
bands. 


14 


4.0 


Figure  II-5.  E^qponentlal  temperature  profile  for  raytrace  in  Figure  II-6. 
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HI.  APPROACH 


Th«  prluzy  motivation  of  this  roaaaroh  effort  was  in  support  of  the  DARXETBS 
submarine  mounted  flir  systesn  being  developed  by  MOSC.  The  system  is  intended 
to  operate  close  to  the  ocean  surface.  Hcerever  for  initial  testing,  it  was 
mounted  or.  a  teat  stand  approximately  30  metres  above  the  ocean  at  NOSC  b'lilding 
IS  on  Point  Lome. 

A.  Characterisation  of  Image  Distortions 

Hm  initial  approach  of  the  research  effort  was  aimed  at  characterising 
image  distortions  obtained  fay  the  OABXEYES  system.  Targets  of  opportunity  (ships 
at  sea  and  the  Coronado  Islands)  wore  observed  during  Santa  Ana  periods  when 
strong  elevated  optical  ducting  oonditione  were  present.  I'ur  this  purpose,  standard 
radiosonde  transmitters  were  modified  to  respond  only  to  tsnpereture  and  pressure. 
The  cepsbility  was  developed  to  routinely  launch  tha  halloons  from  a  boat  at 
saa  along  the  propagation  path.  Standard  surface  meteorological  equipment  wee 
also  installed  at  building  15  to  measure  temperature,  relative  humidity,  pressure, 
visibility  and  aerosol  sise  distributions. 

A  ooiputer  program  dssignated  IMAGE  (Appendix  B)  vaa  developed  to  model 
image  distortion  using  the  measured  temperature  profiles.  Figure  IIZ-1  is  a 
photograph  of  a  disturted  image  of  a  ship  taken  from  a  distanca  of  betwean  10 
and  20  Ion,  during  a  1971  Santa  Ana  oondltlon.  A  temperature  inversion  layer  was 
present  at  a  height  of  approximately  15  metres  ovsr  the  ocean.  In  Figura  111-2 
Is  shown  an  example  of  an  attoupt  to  modal  tlM  distorted  iamge  using  the  IMAGE 
program  and  an  assumed  temperature  gradient  of  .8^C/metre  between  altitudes  of 
15  and  20  mecres.  The  temperature,  T,  and  modified  refractlvlty,  M,  profiles 
are  shown  on  the  left  side  of  the  figure.  (The  modified  refractlvlty  Is  given 
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Figure  III-l.  Distorted  unaige 
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by  M  ■  N  •f(h/a)10^,  yAMX9  h  is  h«lght  above  the  siurface  and  a  la  the  earth's 
radius.  M  is  used  instead  of  N  so  that  the  earth  can  be  represented  by  a  plane 
in  a  raytrace  diagram.  The  fictitious  Increase  vlth  height  of  refractivlty 
provides  that  the  difference  between  the  curvature  of  the  earth  and  the  rays 
remains  the  seine.) 

Shown  on  the  right  side  of  Figure  II1-2  is  a  simulated  ship  configuration 
and  its  image  synthesised  by  raytrace  with  the  indicated  refractivlty  profile. 

The  height  of  the  observer  is  taken  to  be  10  metres,  the  ship's  height  is  20  ii«tres, 
and  the  range  to  the  ship  is  15  km.  A  good  comparison  is  obtained  between  the 
synthesised  ship's  image  and  that  shown  in  Figure  l'IX-1.  The  raytrace  diagram 
used  to  construct  the  distorted  image  is  shown  in  Figure  III-3.  From  a  distance 
of  IS  km,  the  rays  emitted  from  portions  of  the  ship  above  about  3  metres, 
reaching  an  Observer  at  10  metres,  appear  to  be  coming  from  two  different  directions. 
This  results  in  a  two-image  "superior  mirage."  The  upper  image  appears  inverted 
and  elongated  while  the  lower  image  is  i^right  and  slightly  depressed. 

This  describes  the  method  by  which  it  was  initially  Intended  to  characterise 
the  images  obtained  by  the  DABKEXES  system.  However,  during  the  period  (Oct  77- 
Feb  78)  idten  the  system  was  operated,  Santa  Ana  occurrences  were  rare  and  the 
tengwrature  inversions  were  usually  below  the  sensor  height  as  illustrated  in 
Figure  III-4.  For  this  case,  only  the  rays  from  the  lower  portions  of  the  object 
are  bent  by  the  refractivlty  gradient  of  the  inversion  layer.  The  major  portions 
of  the  image  are  therefore  wdlstorted. 

When  elevated  optical  ducts  are  present,  image  distortion  can  occur  only 
when  the  object  and  the  detector  are  in  or  near  the  teirperature  inversion  layer. 

This  is  Illustrated  in  Figure  XII-5  for  the  case  when  the  detector  end  object 
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are  both  slightly  below  the  teroperature  inversion  layer.  Rays  from  the  object 
are  bent  downward  by  the  refractivity  gradient  and  a  second  image  appears  at  the 
detector  inverted  and  at  an  elevated  position. 

During  only  one  mild  Santa  Ana  event  (2  December  1977)  v^s  an  inversion 
layer  observed  near  the  required  30  metres  height.  Figxare  1II-6  is  a  photograph 
taken  during  this  period  showing  a  visually  distorted  image.  The  photograph 
weis  taken  from  a  location  approximately  3  Inn  north  of  the  DARKEYES  facility 
and  from  the  Seuae  height.  The  ship,  determined  to  be  an  LST  at  an  lanknown 
distance,  was  not  detected  by  the  DARKEYES  system.  However,  this  may  have  been 
because  the  ship  was  at  too  large  a  distance  to  have  been  within  the  detection 
sensitivity  of  the  system. 

At  the  end  of  February  1978,  the  DARKEYES  system  was  taken  to  sea  aboard 
a  submarine  for  TECHEVAL,  and  considerable  time  was  spent  in  attempting  to 
obtain  an  alternative  thermal  imaging  system  vhich  could  be  operated  close  to 
the  ocean  surface.  Several  ccminerclally  available  thermal  imaging  systems 
which  use  pj'roelectrlc  vidicon  detectors  were  evaluated.  However,  those  with 
cost  within  the  budget  limitations  of  the  project  were  determined  not  to  have 
sufficient  resolution  and  sensitivity  to  resolve  distant  targets. 

B.  Laser  Transmission  Link 

As  an  alternative  means  of  investigating  ducting  phenomena  at  IR  wave¬ 
lengths  close  to  the  ocean  surface,  the  following  approach  was  adopted.  An  infra¬ 
red  laser  and  detector  were  positioned  such  that  a  line  of  eight  path  existed 
between  then  at  low  tide  (Figure  III-7) .  Here  it  is  assumed  that  at  low  tide, 
the  beam  path  is  sufficiently  elevated  from  the  ocean  surface  so  that  at  no  time 
it  passes  through  regions  of  strong  refractivity.  This  is  probably  a  good  assun^- 
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Figure  III-6.  Photograph  of  ship  taken  in  the  presence  of  a  temperature  inversion 
at  30  m. 


Figure  IIZ-7.  Receiver  In  line  of  slghc  path  of  laser  beam  at  low  tide. 


Figure  I1I-8.  Receiver  In  earth's  shadow  at  high  tide  (T  =  T  ). 

AxJT 
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tlon  for  paths  which  remain  above  a  height  of  about  1  metre,  since  the  major  | 

portion  of  the  air-sea  ter^terature  difference  is  accounted  for  within  the  first  j 

i 

10  cm  beyond  the  interface.  ! 

As  the  tide  rises,  the  sea  siirface  interposes  itself  between  the  source  and  | 

deteccor.  If  no  strong  refractive  effects  exe  present  close  to  the  surface 

(T  ■  T  ) ,  the  beam  would  siitply  be  extinguished  when  the  source  is  eclipsed 
S  A 

(Figure  III-8) .  if  however,  there  is  a  strong  decrease  in  refractivity  with 

height  (T  <  T  ) ,  the  rays  which  pass  close  to  the  surface  would  be  bent  down 
S  A 

around  the  earth  reaching  the  detector  (Figure  III-9) .  Simil2u:ly,  if  there  is  a 
strong  increase  in  refractivity  with  height  (T  >  T  ) ,  the  rays  passing  close  to 
the  surface  would  be  bent  \sp»axd,  causing  loss  of  the  signal  sooner  than  ea^jected 
(Figure  III-IO) . 

Therefore,  by  establishing  an  infrared  transmission  link  close  to  the  ocean 
surface  and  monitoring  signed  strength  as  a  fvinction  of  water  height,  the  refract¬ 
ivity  of  the  atmosphere  just  above  the  air-sea  interface  nay  be  investigated. 
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Figure  III-IO.  Receiver  below  refracted  laser  beam  at  high  tide  (T 


tv.  DESCRIPTION  OP  CO^  lASER  TRANSMISSION  EXPERIMENT 

A.  Transmlsaion  Path 

In  any  experiment  designed  to  measure  small  variations  in  the  refraction 
of  lights  a  long  path  length  is  desirable.  In  addition,  the  propagation  medium 
should  be  as  uniform  as  possible,  with  minimal  effects  due  to  wind  or  wave  motions. 
For  these  reasons,  as  well  as  accessibility,  we  chose  to  conduct  our  measurements 
over  an  8.5  km  path  stretching  from  southeast  to  northwest  over  San  Diego  Bay 
(Fig\ire  IV-1) . 

The  laser  transmitter  was  contained  in  a  portable  van.  It  was  located  on  a 
sand  spit  at  the  San  Diego  Gas  and  Electric  Plant  in  Chula  Vista.  The  laser 
receiver  and  electronics  were  housed  in  an  equipment  hut  situated  on  a  concrete 
pier  at  the  Nave^.  Anphiblous  Base  in  Coronado,  since  the  path  was  out  of  the 
main  channel,  uninterrupted  transmission  could  be  maintained  for  long  periods  of 
time  with  minimum  interference  from  boats. 

B.  Experimental  Apparatus 

1.  General  Description 

The  experimental  apparatus  is  shown  schematically  in  Figure  IV~2.  The 
inl4.ctred  radiation  emitted  by  a  00^  laser  was  chopped  and  directed  toward  the 
receiv^^g  site.  There,  it  was  collected  by  a  cassegralnian  telescope  and  con¬ 
verted  to  an  electrical  signal  by  a  pyroelectric  detector.  The  chopper  reference 
signal,  «^lch  indicates  when  the  chopper  aperture  is  open  or  closed,  was  simul¬ 
taneously  transmitted  by  means  of  a  microwave  radio  link.  This  allowed  a  lock-in 
amplifier  to  be  used  for  synchronous  detection  of  the  laser  signal,  with  a  result¬ 
ant  increase  in  system  sensitivity. 
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Figure  IV-1.  Transmis.iion  path  over  San  Diego  Bay  for  CO.  laser  transmission 
experiment.  ^ 
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2.  LM«r  Tr«naadtt«r 

Tha  main  ccnqponant  of  tha  laaer  tranamittar  was  a  GTE  nodal  941S  CO^ 
lasar.  It  radlatad  3  watts  CW  at  10.6  ym  with  a  long  tam  aaigplltude  stability 
of  batter  than  5%.  Tha  00^  laser  was  rigidly  mounted  on  an  aluminm  plate  along 
with  a  Spactra-Physlcs  modal  132  2  tsM  helium-neon  laser  and  a  Heaver  model  T-16 
16X  gunscopa  (Plgura  lV-3) . 

The  lasers  and  tha  gunscopa  were  successively  sighted  in  at  10  m,  200  m  wd 
1  km.  Small  adjustments  in  their  positions  were  made  until  the  optical  axes  of 
all  three  were  parallel.  This  arrangement  facilitated  aiming  the  CO^  laser  at  the 
receiving  site  since  it  provided  visual  references  for  the  operators  on  both 
ends  of  the  path.  Hhan  the  r«>.eivlng  site  was  placed  within  tha  crosshairs  of 
the  gunscopei  the  hellum-neoa  lasar  was  easily  visible  to  an  observer  at  that 
location  (even  in  the  daytime) . 

In  order  to  achieve  good  beam  pointing  stability,  the  laser  plate  was  mounted 
across  3  inch  I-beams  which  were  bolted  directly  to  the  frame  of  tha  van.  Tha 
van  itself  was  raised  on  concrete  blocks.  Mo  appreciable  pointing  jitter  was 
noted. 

The  Infrared  laser  beam  was  chopped  at  40  Hz  by  a  rotating  blade  chopper 
salvaged  from  a  surplus  Laser  Precision  model  FKP-345  radiometer  probe. 

3.  Laser  Receiver 

The  chopped  infrared  radiation  was  collected  by  a  Molectron  model  PT 
cassegralnian  telescope.  It  had  an  effective  focal  length  of  200  nan,  an  effective 
aperture  number  of  f/2.8,  and  a  full  angular  field  of  view  of  1.42  degrees.  The 
infrared  radiation  was  detected  by  a  Molectron  model  P4-45  pyroelectric  detector. 
It  was  fitted  with  an  anti-reflection  coated  germanium  window  in  order  to  reject 
background  light  outside  the  8-14  ym  band.  The  responsivity  of  the  detector 

4 

was  on  the  order  of  10  volts/watt. 
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Th*  •IttCtrical  signal  from  the  detector  ms  monitored  by  a  Princeton  Applied 
Pesearoh  model  5204  lock-ln  analyzer  In  the  vector  voltmeter  &iode>  The  magni¬ 
tude  of  signal  ms  observed  on  a  front  panel  meter  and  continuously  recorded  on 
a  Mosely  model  680  strip  chart  recorder. 

The  telescope-detector  combination  ms  mounted  on  an  adjustable  tilt-plate 
vhich  could  be  located  on  any  desired  support  at  the  receiving  site.  In  Figure 
IV- 4,  It  Is  shorn  In  place  on  a  specially  constructed  platform  which  suspended 
the  assembly  close  to  the  ocean  surface. 

4.  Microwave  Reference  link 

In  order  to  employ  synchronous  detection.  It  was  necessary  to  provide 
a  chopper  reference  signal  to  the  lock-ln  amplifier.  This  was  done  by  means 
of  the  microwave  reference  link  shown  In  Figure  IV-5. 

A  2.9  GRz  Fairchild  phase- locked  oscillator  ms  amplified  to  about  0.5  watt 
by  a  Hewlett-Packard  model  491C  traveling  mve  amplifier.  The  THA  was  modulated 
by  the  chopper  reference  signal  and  fed  to  a  20  dB  gain  microwave  horn  for  trans¬ 
mission. 

The  signal  was  received  by  a  similar  horn,  amplified  by  an  EH  Labs  model 
2900  preamplifier,  and  then  fed  to  a  Mlcrotel  model  1200  frequency  tracking 
receiver.  The  chopper  signal  ms  extracted  from  the  carrier  with  a  Hewlett- 
Packard  model  423A  diode  detector  and  processed  with  a  comparator  amplifier 
in  order  to  produce  a  clean  reference  signal  for  the  lock-ln  amplifier. 

C.  Meteorological  Measurements 

During  the  course  of  a  run,  the  following  met  measurements  were  made 
at  the  receiver  site.  Air  and  sea  tesperatures  were  measured  to  within 
using  themiometers.  The  air  tesperature  was  generally  measured  at  about  2  m  above 


35 


Figure  IV“4,  Phol^ogreph  of  lelescope-detector  conblna'tion  suspended  close  to 
the  ocean  surface. 
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Figure  IV— 5.  Schematic  diagram  of  microwave  reference  link. 


the  ocean  s\ir£ace.  Wind  speed  was  measured  to  within  aJjout  nph  using  a  hand 
held  instrument.  The  height  of  the  detector  above  the  water  line  was  measured 
to  within  about  j^5  cm  using  a  dipstick. 

D.  Measurement  Procedure 

Each  time  data  were  taken  i  measurements  were  made  in  the  following  manner. 
The  microwave  horns  were  set  up  and  ooaorsely  pointed  at  the  opposite  site.  The 
telemetry  was  then  turned  on  and  reception  of  the  reference  signal  verified. 

The  helium-neon  laser  was  then  turned  on.  Zn  general,  it  was  visible  at 
some  height  at  the  receiving  site  which  was  dependent  on  the  meteorological  con¬ 
ditions  and  the  tide.  Realignment  of  the  laser  plate  was  usually  not  atten^ted 
unless  the  beam  was  not  visible  at  all. 

An  eyepiece  was  substituted  for  the  pyroelectric  detector  in  order  to  visu¬ 
ally  align  the  collector  telescope.  The  detector  was  replaced  and  the  CO^  laser 
turned  on.  Since  the  beam  divergence  of  the  helium-neon  laser  was  1  mr  and  that 
of  the  00^  laser  was  4  mr,  the  infrared  signal  was  almost  always  observed  after 
this  process.  The  lock-in  amplifier  was  run  with  a  time  constant  of  100  seconds 
in  order  to  smooth  out  signal  flucttiatlons .  Meteorological  measurements  were 
made  at  the  receiving  site  approximately  every  half  hour  during  the  run. 
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DATA 


I 

A.  General  Behavior  | 

i 

Over  the  course  of  three  months ^  data  was  taken  on  eighteen  separate  | 

days.  During  all  days  but  the  first,  the  detector  had  to  be  located  significantly  \ 

I 

higher  thw  would  be  e^^pected  for  line  of  sight  contact  in  the  standard  atmosphere.  | 

Throughout  these  runs,  the  sea  tesperature  generally  exceeded  the  air  tesperature 
by  several  degrees  centigrade. 

For  data  taken  during  a  rising  tide  (Figure  V-1) ,  the  laser  signal  (while 
it  was  detectable)  exhibited  no  shaurp  fades  or  fluctuations.  As  the  tide  rose 
and  contact  was  lost,  the  signal  dropped  several  orders  of  magnitude  In  just  a 
few  minutes.  When  this  occurred,  the  detector  was  raised  ip  until  the  signal 
was  again  acquired.  The  procedure  was  repeated  until  high  tide  was  reached. 

Data  taken  during  falling  tides  (Figure  V-2)  vnis  similar  except  that  the  detector 
was  successively  lowered  each  time  the  signal  emerged  from  the  background. 

B.  Analysis 

Let  an  optical  source  at  height  h  and  a  receiver  at  a  different  height 

8 

h^  be  separated  by  a  distance  R.  Further,  as  shown  in  Figure  V-3,  let  the 
source  and  receiver  be  positioned  so  that  their  optical  horizon  distemces  d^  and  d^ 
respectively,  coincide.  Therefore  d  ■<-  d  ■=  R  and  at  this  point,  the  source 

8  3- 

is  just  detectable  by  the  receiver.  The  distance  to  the  optical  horizon,  d  , 

o 

as  seen  by  an  observer  at  height  h  through  a  refractive  index  gradient  dn/dh 

o 

is  given  by  (Reed  and  Russ'^ll,  1964) 

<2o 

K  l+a(dK/dh) 

where  a  is  the  radius  of  the  earth. 
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v-3.  Diagram  of  horizon  distances  for  source  and  receiver  at  heights 
h^  and  h^  respectively. 


Then  from  (1) 


d  *  k  \/h~ 

S  Vs 

(2a) 

and 

d  ■  k 
r  V  r 

(2b) 

where 

k  «  J  2a 

(3) 

If  l+a(dn/dh) 

or 

d  +  d  =  R  -  k  (Vh~  -kST) 
s  r  ^  s  V  r 

(4) 

If  R  is  fixed,  h  and  h  may  be  e^^ixnentally  determined  from  the  water 

8  f 

height  at  the  time  the  observer  just  loses  sight  of  the  source.  Then  the  horizon 

distance  for  an  observer  at  any  height  h  (for  the  same  llne£u:  refractive  index 

o 

gradient  as  for  which  and  h^  vexe  measured)  is  given  by 


«here  the  units  of  d  are  those  of  R. 

o 

The  values  of  h^  and  h^  measured  during  14  different  eclipses  or  appearances 

of  the  signal,  along  with  the  corresponding  values  of  k  zure  shown  in  Table  V-1. 

Also  shown  are  the  oorrespondlng  values  of  measured  sea  and  air  teiqperature  (air 

ten^rature  measured  at  fixed  height  on  pier).  With  the  exception  of  3  events, 

the  air  teit^rature  was  always  less  than  the  sea  temperature  (these  are  end-point 

measurements  and  are  probably  not  representative  of  average  values  over  the  path) . 

The  average  value  of  k  for  the  14  events  is  2.15  Then  from  (3)  the 

“7  -1 

refractive  index  gradient  can  be  determined  to  be  2.76x10  m  .  From  equation 
(III-3)  the  ten^rature  gradient  corresponding  to  this  gradient  of  refractive  index 
is 


I?  ^-1 


where  N  is  given  by  III-2.  The  average  value  of  T  for  the  data  set  is  22.7  C. 
For  the  measured  values  3T/B  z  % -.3®c/roetre. 
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Table  V-1 


Date 

T 

T  4 

rc) 

hg 

(m) 

(m> 

k 

-■k 

knric 

11  Jul 

21.7 

19.9 

3.66 

4.89 

2.07 

11  Jul 

22.1 

20.0 

3.24 

5.13 

2.10 

12  Jul 

22.0 

21.7 

3.47 

4.70 

2.12 

13  Jul 

23.2 

24.8 

3.09 

4.32 

2.23 

17  Jul 

23.3 

18.8 

3.55 

4.78 

2.10 

20  Jul 

22.7 

21.9 

3.26 

5.15 

2.10 

20  Jul 

22.6 

21.9 

3.55 

4.78 

2.10 

21  JUl 

22.7 

20.2 

3.02 

4.91 

2.16 

24  Jul 

22.2 

22.3 

3.02 

4.24 

2.25 

27  Jul 

22.9 

21.5 

3.24 

4.47 

2.18 

28  Jul 

23.2 

21.9 

3.24 

4.47 

2.18 

4  Aug 

23.6 

23.8 

3.40 

4.63 

2.14 

23  Aug 

23.2 

21.5 

2.94 

4.83 

2.19 

23  Aug 

23.1 

21.8 

3.17 

4.40 

2.20 
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Tlj«  distance  to  the  geometrical  horizon  is  found  from  (1)  for  dn/dh«0  to 
be  d  »  S.Siyh 

•.rm.  g 

where  h  is  in  metres  and  d  is  in  kilometres.  The  results  of  this  study  show  that 
for  negative  temperature  gradients  of  a  few  tenths  of  a  degree  centrigrade,  the 
distance  to  the  infrared  (or  optical)  horizon  (d^%  2.15^h)  is  only  about  60% 
of  tlie  distance  to  the  geometrical  horizon. 

The  distance  to  the  optical  horizon  for  a  wide  range  of  tenperature  and  index 
gradients  is  shown  in  Figure  V-4.  It  is  evident  that  the  horizon  shortens  relative¬ 
ly  slowly  for  temperature  gradients  more  negative  than  that  of  the  standard 
atmosphere.  For  positive  gradients  however,  the  corresponding  increase  in  horizon 
distance  is  much  more  rapid.  This  indicates  that  temperature  gradients  of  similar 
magnitude  but  opposite  sign  to  those  observed  in  this  experiment  would  result  in 
greatly  extended  horizons. 

Listed  in  Appendix  C  is  a  conputer  program  called  BAYFRDP.  It  may  be  used 
to  compute  the  various  parameters  related  to  a  source  or  observer's  horizon  for 
a  constant  gradient  of  index  of  refraction. 
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Figure  v— 4,  Distance  to  the  horison  va»  iJsd^A,  gradient  and  teoiperature  gradient 
for  an  observer  at  a  height  of  1  m. 


VI.  CONCLUSIONS  MiO  BECO^a^ENIIAXIONS 


Because  o£  the  prevailing  meteorological  conditions  which  existed  In  San 
Diego  Bay  during  the  measurement  period,  this  study  was  limited  to  data  acquired 
\dien  the  bay  water  was  wanner  than  the  air.  For  air  temperatvure  gradients  less 
than  -l^C/m,  it  was  demonstrated  that  a  severe  shortening  of  the  optical  and 
infrared  Iv^xizons  consistently  occur.  Also,  calculations  indicate  that  if  gra~ 
dients  of  6ii..ilar  magnitude  but  of  opposite  sign  occur  close  to  the  ocean  sur¬ 
face,  greatly  extended  horizons  are  to  be  expected.  However,  most  of  the  advan¬ 
tages  gained  by  extended  horizons  cotild  be  nullified  by  the  effects  of  absorptior. 
due  to  water  vapor,  in  any  case,  since  in  many  of  the  Navy's  operational  areas 
(especially  during  winter]  the  sea  is  wanner  than  the  air,  the  detrimental  effects 
of  shortened  horizons  and  distorted  images  are  more  likely  to  be  encountered. 

In  addition  to  the  above  effects  to  be  «qpected  on  passive  surveillance 
devices,  laser  systems  (target  designators  aivi  trackers)  operated  close  to  the 
ocean  surface  could  be  subject  to  large  errors  in  pointing  and  tracking.  This 
could  result  from  strong  refractivity  gradients  which  exist  near  the  air-sea  Inter¬ 
face. 

The  quantitative  studies  made  here  have  been  limited  to  infrared  wavelengths 
propagating  close  to  the  ocean  surface.  However,  the  visual  evidence  presented 
indicates  that  elevated  ducts  could  also  cause  severe  image  distortions  for 
imaging  systems  operated  at  other  heights  (aircraft,  ship  mast,  etc.) . 

Based  on  the  results  of  this  study,  it  is  resonnended  that  a  climatological 
data  base  of  air-sea  teroperattire  differences  for  operational  areas  be  asseinbled. 
This  would  aid  shipboard  oommanders  In  assessing  the  effects  of  refractivity 
gradients  on  electro-optical  systems'  performance.  For  exasfile,  a  library  could 
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be  developed  to  ascertain  the  distance  to  the  horizon  for  a  particular  shipboard 
electro-optical  syatem  at  any  location  and  tine  of  year. 

Peal  tine*  shipboard  optical  refractivity  sensors  also  need  be  developed  to 
enhance  predictive  capabilities.  The  conputer  programs  developed  in  this  investi¬ 
gation  provide  a  means  of  predicting  image  distortion  in  the  presence  of  refract¬ 
ivity  gradients.  These  programs,  in  conjunction  with  real  time  meteorological 
measurements,  «ould  enable  simulated  Images  of  expected  targets  to  be  displayed. 
This  would  aid  electro-optical  systems  operators  in  identl^ing  targets  under 
degraded  optical  conditions. 


REFERENCES 


1.  FLEAGLE,  R.  G.,  "The  Optical  Measurement  of  Lapse  Rate",  Bull.  Amer.  Met. 
Soc.,  31,  51  (1950) 

2.  FLEAGLE,  R.  G.,  "The  Temperature  Distribution  Near  a  Cold  Surface",  Jour, 
of  Met.,  13.  IGO  (1956) 

3.  FLEAGLE,  R.  G.  and  BUSINGER,  J.  A.,  "An  Introduction  to  Atmospheric 
Physics",  Academic  Press,  CHAPTER  VII,  (1963) 

4.  KERR;  D.  E.,  "Propagation  of  Short  Radio  Waves",  McGraw-Hill  (1951) 

5.  REED,  H.  R.  and  RUSSELL,  C.  M.",  "Ultra  High  Frequency  Propagation", 
Boston  Technical  Publishers,  Chapter  3,  (1964) 

6.  SELBY,  J.  E.  A.,  SHETTLE,  E.  P.  and  MC  CLATCHEY,  R.  A.,  "Atmospheric 
Transmittance  from  .25  to  28.5  urn:  Supplement  Lowtran  3B  (1976)", 

AFGL,  Paper  No.  AFGL-TR-76-0258 


49 


APPENDIX  A 

OOMPUIER  PROGRAM  RAYTRACE 


This  appendix  describes  a  FORTEAN  computer  program  that  calculates  and  dis¬ 
plays  ray  trajectories  for  arbitrary  but  piecewise  linear  modified  refractlvity 
versus  height  profiles. 

Program  PAYTEACE  is  a  FOETRAH  IV  progrtun  designed  for  use  on  the  Data  GenereQ. 
Mova  SCO  mlnlconputer  system  using  the  operating  system  EDOS.  Both  input  and 
output  of  the  program  is  via  the  T^tronix  4012  interactive  graphics  display 
terminal.  The  purpose  of  the  program  is  to  calculate  and  display  ray  trajectories 
on  a  height-versus-range  plot  for  axi  arbitrary  but  piecewise  linear  refractivity- 
versus-height  profile.  Other  inputs  to  the  program  are:  maximum  range,  range 
irtcrement,  maximum  height,  source  height,  initial  launch  angle,  final  launch 
angle,  and  number  of  rays.  The  program  uses  Snell's  law  with  small  angle  approx¬ 
imations  to  compute  the  ray  trajectories  based  on  an  input  modified  refractlvity 
profile  specified  in  M-units.  Modified  refractlvity  in  M-units  is  related  to 
refractlvity  in  N-units  by  the  relation 
K  «  N  +  0.048  h 

where  h  is  in  feet.  This  appendix  will  not  attempt  to  e:q>lain  the  inner 
workings  of  the  program  or  the  formulations  employed,  but  rather  will  outline 
the  operation  of  the  program. 

Figure  A-1  shows  a  sample  of  the  proopt  and  input  sequence  for  entering  the 
modified  refractlvity  profile  that  occurs  at  the  beginning  of  the  program.  Note 
that  the  required  levels  must  be  entered  in  order  of  descending  height  and  must 
end  at  the  zero  height  level.  A  nuuclmum  of  50  levels  may  be  entered.  After  the 
zero  height  level  has  been  entered,  the  program  clears  the  screen,  plots  the 
EP-unlts  versus  height,  and  prompts  for  the  remaining  inputs  as  illustrated  by  the 
example  of  Figure  A-2.  All  of  the  remaining  inputs  are  self  explanatory  except. 


INPUT  PROFILE  FRON  HIGHEST  TO  LOUEST  UZTH  THE  FORM 
HCZQHTCFEET)*  N  UNITS 


perhaps,  l^e  range  increment  vdilch  should  be  selected  based  on  the  maximum  rzuige 
and  the  ccmiplexlty  of  the  M-unlt  profile  to  give  smoothly  curving  ray  trajectories 
Under  most  conditions  a  range  Increment  of  one-twentieth  the  maximum  range  Is  a 
good  choice. 

Figure  A-3  shows  the  final  ray  trace  plot  for  the  inputs  shown  in  Flg\ires  A-1 
and  A-2.  After  each  ray  trace  plot  is  generated,  an  entry  of  "0"  (for  old  Profile 
followed  by  a  carriage  return  will  direct  program  flow  to  the  entries  illustrated 
in  Figure  A-2.  An  entry  of  "N"  (for  New  Profile)  followed  by  a  carriage  return, 
will  direct  program  flow  to  the  M-unit  profile  inputs  illustrated  in  Flgvire  A-1. 

A  complete  FOFXRAN  program  listing  of  program  PAYTPACE  and  the  one  subroutine 
TRACE  are  included  here  following  the  figures. 


ce  Plot  for  the  Inputs  Specified  in  Figures  A-1 


prograA  raytrace 

X<50^,R«(SG>.&HLIP(6)tHLlH<7),L&L(6>tLeA(13> 
COPMOM/B/  QA(5C5 ,YAC50> ,NPAX 
DATA  OWLIW/'C.  ,l5r!.,200.t5CC..10CO,.200C./, 

HL1K./5CC.  .iGC'c..2cri..5grg.,ioooc.,2CLCo.  ,40pca./ , 

'  LEA/A,L«T,I,T,UfD»E»  »F»EtEt 
CALL  1N1TT(460) 


2  0 

FULL 
COFIPLETE 


30 

35 


4C 


60  TO  35 


OWLIPCI)  ;T0TAL  M 

II  »  K«mN-APi00<RM'*IN,0K/5.)  ;''1M»'UW 
*H(?) 


1 
2 

3  LEA/'A'.'L' 

CALL  1N1TT(460) 

,  TYPE  ‘‘input  profile  FRO*' 
1”HEIRHT<FEET>t  «  UNITS 
00  lb  1=3, 5C 
ACCEPT  ••  »  ".HdKRWd) 
IF<H(I) .tO.C.O)  GC  TO  ‘ 

1C  CONTINUE 

TYPE  "  PROFILE  ARRAY 
60  TO  5 

PROFILE  READ  IN 
2Q.J4N/g(«l  - 

RPMIN  s  RN(NMAX) 

RPNAX  =  RNKIN 
00  25  I=3,NKAX 

Rjei  s  R«(i)  .  - 

IF(RPI.LT.RPMIN)  RNNlNsRKI 

if<rmi.6T.ppmax5  rnmax=rni 
25  CONTINUE 

OELH  =  RMKAX-RP.MIN 
00  3C  1»1,5 
iFCDELK.LT.OMLIPd)) 

CONTINUE 
1*5 
OP 

PPPTN  -  n. 

HPAX*H(?) 

00  40  1*1,7 

IF<HPAX,LE.HLIP<1))  60  TO  45 
CONTINUE 
1*7 

45  HNAX  =  HLIP^Cn 
C  DRAW  VERTICAL  AXIS 
47  CALL  NEWPAG 
IY»2C0 

CALL  P0VAeS-(190,1C0X 
CALL  ORWABS  (20C,irO) 

DO  50  1*1,5 

CALL  DRWABS  (20C,1Y) 

CALL  DRWABS (19C,1Y) 

CALL  POVABS  (200,1Y) 

50  1Y*1Y+1PC 
C  DRAW  HORIZONTAL  AXIS 

CALL  POVABS  UCC, 90 
CALL  DRWABS  (2CC, 1CC) 

IX*26C 
DO  60  1*1,5 
CALI.  OJlMABACI.X,iaCU 
CALL  ORWAPS  (IX  ,90) 

CALL  *0VAPS <IX  ,10C) 

60  IX*IX-»6? 

C  PLOT  P  VERSUS  HEIGHT 

UPM*300./DW  ;SCR££N  UNITS 
UPH  =  5Ca./h»’AX  ;SCFEEN  UMTS 
IX  =  2D0-»  (  R  •'C?  )-FPNIN)  *UPP 
lY  *  100*H(7)*UPH 
CALL  MOVABS  (IX  ,1 Y) 

00  7C  I*4,NPAX 
IX  *  20D'»(RP(I  )«F»f''^IN)*UP'' 
lY  »  10r+H(I)*UPh 
70  CALL  DRWAFS  (JX  ,IY) 

1Y»41C 
00  75  1*1  , < 

CALL  POVAPS  (f 0  ,IY) 

CALL  ANPODE 


-  HI6HES>'  TO 
<15>",”<15>** 


LOWEST  WITH 


'T'/ 

THE  F0RM<15>'*, 


RE-ENTER  PROFILE*' 

SCALE  PROFILE  FOR  PLOTTING 


CHANGE  ON  AXIS 
P  ON  AXIS 


PER 

PtC 


P  UNIT 
FOOT 


A6 


KmTE<10»10C)  LBL(I>  . 

IY«lV-20 

F0RMAT(1X,A1) 

lt*590  . 

IHLABEL  »  HFAX 
IHINC  «  H^^AX/5 
DO  80  I*1,6 

CALL  KOVAB.S..112n.t-IJL3L _ 

CALL  ANMODE 
WRITECIQ.IZC)  IHLABEL 
IY*IY-100 

1HLABE.L- =  _l.HLAB.ELililliiljC. 
FORMATdX.lS) 

CALL  FIOVAeS(17G,60> 

CALL  ANMODE 

.J!Rlli*PM»LIJil _ 

HHAX>PNM1N4DM 
WR1TE(10,UC)  MMIH 
FORMAT (IXjlA) 

-CAL4.  MOVABSLA70*Aai. _ 

CALL  ANi*ODE 
WRITEdOfUC)  MMAX 
CALL  MOVABS (200*20) 

-CALL  A.liMXU)£ - — 


WRITE  (10*160 
160  FORMAT  (7K«’*R  UNITS'* 
1AY  TRACE  PLCT'y/) 


//////'  FOLLOWING  PARAMETERS 


'  THE  R 


1  DELR*"  MAX  HEIGHT  (FEET)  ■  “.HRMAX, 

2  *•  SOURCfe  height  (FEET)  *  ",HSOURCE,**  INITIAL  LAUNCH  ANGLE 
3=  ".ALPHAS,"  FINAL  LAUNCH  ANGLE  (MRAD)  =  ".ALPHAT,"  NUMBER 

■■■  4-  »-*,NRAYS _ ^ . . . . . . . 

SCALE  PROFILE  FOR  RAYTRACE  PLOT 
$X  «  800/RMAX 
SY-500/HRMAX 

- - - - - 

YA(2)  •  500. 

YA(3)  ■  H(3)*SY 


QACI)  •  3.6E*8«SY0SX2 
00  210  I>A*NMAX 
YA(1)  •  H(X)*SY 

21C  OA(I)  «  (RM(I-1)-RM(I))/(H(I-1)-H(I))*SY0SX2*1  .E-G 

SY-08X  »-S-Y-/-&X- . - . . - . - . - . .  - 

AS  «  ALPhAS*SY0SX*1 .E-3 
AT  «  ALPHAT<»SY0SX*1  .E-3 
YS  «  HS0URCE*SY 
DELA  *  (AS“AT)y(NRAYS-1) 

AL  «  AS 


(MRAD) 
OF  RAYS 


AT  «  ALPHAT<»SY0SX*1  .E-3 

YS  «  HS0URCE*SY 

DELA  *  (AS“AT)y(NRAYS-11 

AL  «  AS 

OR*DELR*SX 

DO  215  I«1,NMAX 

J>NMAX-I 

IF(YS.LT.YA  (J)  )  60  TO  21 
215  CONTINUE 
217  KS»J*1 

C  DRAW  VERTICAL  RAYTRACE  AXIS 

CALL  NEWPAG . . 

CALL  MOVABS  (UO,  ICO) 
CAL^^DRWABS  (150*100) 

DO  22C  1=1,5 

CALL  DRWABS  (ISC,  lY) 

CALL  DRWABS  (14C, IV) 
call  MOVABSdSC.IY) 

220  IY*lY-dnO 
ir®48C 

DO  230  1  =  1  ,  13 
CALL  MOVABS  (30  ,IY) 

CALL  ANMODE 
WRITEdO*31C)  LBA(l) 

230  IY«1V-2C 


A7 


510^  f<MMIAT<tX,AU 
iy*590 

lHLABEL«HRn«X 

1HINC>bHRMAX/5. 

90  ?40  1«1,« 

CALL  WOVABS  (70  ,I Y) 

CALL  AN**ODE 
WR1TE{10,12C)  IHLAPEL 
IY»IY-10&  - 

240  IHLA8EL  «  X HLABE L-I H IN C 
C  DRAW  HORIZONTAL  RAYTRACE  AXIS 
CALL  M0VA8S  <150,50 
CALL  DRMABS  CISC,  ICO) 
IX»31C 

DO  25D  I»1,5 
CALL  DRWABS  (IX  ,1  CO 
CALL  DRWABSdX,  90 
CALL  HOVABS  (IX  ,1  CO 
250  IX  »  IX+16G 

CALL  M0VA6S(38C, 

CALL  ANMODE 
HRITE(10,?2C) 

320  FORHATdX, ’range 
IR»C 

-  ■  .  -  ^4A^C*A*AX/-5* . 

IX«120 


20) 


IN  NAUTICAL  NILES') 


60 


DO  260  I«1,6 
CALL  HOVABSCIX, 

CALL  ANNODE 
WRITE  (40*330  IR 
IR  *  IR+IRINC 
IX«IX-»16C 
FORNAT (lx  ,14) 

TRACE  ALL  RAYS 

DO  270  isl.NRAYS 
CALL  TR  A C E  <  YS  ,  C  AL-,0  R ,  K  S ) 
A'L  *  AL^OELA 
CALL  nQVA8S(9CC,0 
CALL  ANHODE 
RfAD(11,5C0)  lOPT 
F0R"'AT(A1) 

CALL  KEVPAG 
IF(10PT*E0. 'N  ')  GO 
IFdOPT.EQ. '0  ')  GO 
CALL  FINITT  (0,830 
END 

SUBROUTINE  TRA  CF  (Y$  ,  ,  •«  9  ,  u 
C0MH0N/B/QA(5C'),  YA(5C)  ,NMAX 
IX«XS 
IY*YS 


260 

330 


570 


5  CO 


TO 

TO 


5 

47 


XS,AS,DX,KS} 


CALL  KOVAPS (IX+1 5C, I Y+1C:) 

K»KS 
A  s  AS 
X  =  XS 
Y  =  YS 
Q  «  Q  A  ( K  } 

1F(AS.GT,C.)  GO  TO  10 
YLIM  *  YA(K) 

60  TO  20 

1C  YLIN  =  YA(K-I) 

20  DELX  s  DX 

3(5  B  *  A-^DELX^GAOCBC. 

IF  (A.EQ.C.)  GO  TO  115 
IF(SIGN(1  .  ,P)  .  NE  .SIGNd  ..,A)  )  GO  TO  40 
DELY  »  7O40.*(e'*-A)*OELX 
60  TO  50 
40  B  «  0. 

OELX  =  -A/Q/60eu. 

OELY  =  3C40,*A*DELX 
1F(A.GT.C.)  GO  TC  42 


A8 


IF(Y+DELY.LT.YLI«)  GO  TO  52 
60  TO  70 

>2  IF (Y+OELf .6T.YLIW)  GO  TO  102 

60  TO  7-0  . 

10  IF(A.6T.O.)  GO  TO  100 

lF<Y+DELY.6EtYLIM)  60  TO  7C 
iZ  S6N  *  -1. 

.  DELY  =  YLLMrY. .  „ 

IFCK.EO.Kt^AX)  GO  TO  6C 
K=K  +  1 

YLIM  =  YA(K) 

,n  p=CCM*gOBT  itZ . *& * DlELXIJL- 

DELX  »  (B-A)/Q/6DPC. 

Q  «  QA(K> 

'0  X  *  X^OELX 

Y  «  Y+OELY 

If  <X  .GE.SOO.^)  GOTO  114 
IF(Y.6E.50C.)  60  T0^120 
IF(Y.LE.O.)  GO  TO  90 
A  »  B 

>0  ix»x 

CALL  DRWABS  (1X+150,1Y*100) 

60  TO  20 
50  A  =  -B 

YLl*'  =  YA  (NfAX  -1  ) 

60  TO  80 

:0  IF  (Y^OELV »LE.YL1H)  GO  TO  70 

32  S6W  »  1- . 

DELY  =  YLIM-Y 
K  =  K~1 

YL1M*YA<X-1 ) 

60  TO  60 

to  lF<X.EQ«cCC.)  GO  TO  120 
X  *  X-DELX 

Y  *  Y-OELY 

DELX  =  Jf&O.-X  . 

GO  TO  30 

15  DELY  =  <CAC.*CB+A)*OELX 
IFCO.LT.C.)  60  TO  117 
YLIie»YAU-1) 

GO  TO  42 
17  YLIM=YA<K) 

60  TO  *1 
20  IX^X 

CALL  ORWAPS  (IX^I 50, IY  +  1CG) 

RETURN 

END 


APPENDIX  B 

COMPUTER  PPOGRAM  IMAGE 


This  appendix  describes  a  FORTRAN  computer  program  that  calculates  and  dis¬ 
plays  optical  image  distortions  created  by  arbitrary  but  piecewise  linear  temper- 
ature-versus-helght  profiles. 

Program  IMAGE  Is  a  fortran  IV  program  designed  for  use  on  the  Data  General 
Nova  800  minicomputer  system  using  the  operating  system  RDOS.  Both  Input  and 
output  of  the  program  is  via  a  Tektronix  4012  interactive  graphics  display 
terminal.  The  purpose  of  the  program  is  to  calculate  and  display  the  optical 
image  that  an  object  would  exhibit  to  an  observer  at  an  arbitrary  remge  and 
height  above  the  earth's  surface  for  arbitrary  but  piecewise  linear  teaqperature- 
versus-height  profiles  of  up  to  four  linear  segments.  A  fixed  object  corresponding 
to  a  free-foxm  ship's  profile  20  metres  high  and  approximately  108  metres 
long  is  used  in  the  program.  Provisions  are  made  in  the  program  for  changes  in 
the  xefractlvity  based  on  the  optical  wavelength  of  concern.  The  image  is  cal¬ 
culated  by  tracing  a  series  of  rays  through  the  atmosphere  for  uniformly  varying 
elevation  angles  at  the  observation  point.  The  one  metre  high  horizontal  "slice" 
of  tlie  object's  form  closest  to  the  height  computed  by  the  ray  trace  at  the  obser¬ 
vation  range  is  then  plotted  at  the  apparent  height  based  on  the  elevation  angle 
at  the  observation  point  assuming  no  refraction. 

Figures  B-1  and  B-2  are  sanple  outputs  from  program  IMAGE.  Figure  B-1  is 
for  no  change  of  temperature  with  height  and  shows  an  image  that  corresponds 
exactly  with  the  object.  Figure  B-2  is  for  no  increase  of  temperature  to 
28  metres,  then  a  rapid  increase  over  the  next  4  metres,  and  no  increase  above 
that  point.  The  image  in  this  case  is  the  undistorted  object  below  and  an  in¬ 
verted  and  slightly  extended  object  at  higher  elevation  angles  created  by  the 
"trapping"  or  bending  downward  effect  of  the  rapid  temperature  increase  in  the 


Progran  Output  for  the  Case  of  a  Strong  Temperatiure  Increase 
:6  to  30  metres. 


26  tr^  30  Mtr«  Interval.  In  Isoth  cases  an  optical  wavelength  of  .50  microns, 
an  observing  height  of  25  metres,  and  a  range  to  the  object  of  20  hn  are  assumed 
Both  tenperatTire  in  degrees  Celsius  and  modified  refractlvity  In  M-units  are 
plotted  and  numerically  listed  versus  height  on  the  output  display.  Figures 
B-3  und  B-4  eure  ray  trace  diagrams  that  show  the  trajectories  in  space  of  the 
rays  that  create  the  images  in  Figures  B**l  and  B-2  respectively. 

Figure  B-5  is  a  sairple  input  session  for  the  inputs  corresponding  to 
Figure  B-2.  First  the  program  prompts  "ENTER  HEIGHT  (H) ,  gSMP  (C) ,  EESCENDIHG 
TO  ZEBO-5  LEVELS  MAX"  after  ^Ich  the  operator  enters  vip  to  5  sets  of  heights 
and  temperatures.  The  last  height  level  must  be  sero.  Next,  the  program 
prompts  "OPTICAL  WAVELENGTH  (KICRCXtS)  -  ",  "OBSERVING  HEIGHT  (M)  -  ”,  and  "RANGE 
TO  OBJECT  (KM)  «  "  after  which  the  operator  inputs  the  required  numeric  values. 
This  concludes  the  input  session  and  the  program  begins  calculations. 

The  program  calculates  the  optical  refractlvity  N  from  the  relation  given 
in  reference  1: 

K  77.6  P  .  0.584  P 

r  Tx" 

>diere  P  is  atmospheric  pressure  in  mb,  T  is  temperature  in  Kelvins,  and  X  is 
the  optical  wavelength  in  microns.  The  program  assumes  P  decreases  with  height 
at  a  rate  of  .12  ob  per  metre  starting  with  a  value  of  1000  at  height  zero. 

The  modified  refractlvity  M,  which  accounts  for  earth  curvature  is  calculated 
according  to  the  relation 
M  »  K  +  0.157  H 
where  H  is  height  in  metres. 

The  ray  trace  equations  used  in  program  IMAGE  utilize  small  angle  approxima 
tione  since  the  angles  of  concern  are  generally  less  than  1  to  2  degrees.  The 


B4 


Figure  B-3.  Ray  Trace  Diagram  Showing  the  Ray  Trajectories  for  the  Image 
Output  cf  Figure  B-1. 


hanqk  zn  KZLonmp» 


Figure  B-4.  Ray  Trace  Diagreun  Shewing  the  Ray  Trajectories  for  the  Image 
Output  of  Figure  B-2. 
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Figure  B-5.  San^le  Input  Session  for  the  Output  Shown  in  Figure  B-2. 


most  iioportant  program  variables,  which  correspond  to  the  program  listing  included 
in  this  appendix  after  the  figures,  are  defined  below: 

H  (1) :  Height  in  metres  descending  with  increasing  I 
T  (I) :  Temperature  in  degrees  Celsius  corresponding  to  H  (I) 

BM  (II  ;  M  values  corresponding  to  H  (I) 

RS  (I)  :  Gradient  xlO®  of  height  to  M  value  for  the  lineau:  segment  below  K  (I) 


IMAX: 

W: 

HO; 

BKM: 

R: 

AS: 

ElBIiA: 

AO: 

Al: 


(I)  ^  10®  H  (I)  -  H  (I  t  1) 

RM  (I)  -  RM  (1+1) 

Total  number  of  levels.  Note  H  (IMAX)  =  0 

Optical  wavelength  in  microns 

Observing  height  in  metres 

Pcmge  to  object  in  Jon 

Range  to  object  in  metres 

Elevation  angle  at  observing  point  for  each  ray  in  radians 
Incremental  value  of  AS  between  successive  rays 
Elevation  angle  at  beginning  of  each  step  along  ray  in  radians 
Elevation  angle  at  end  of  each  step  along  ray  in  radians 


RTOT:  Range  along  ray  in  metres 


RMO;  M  value  at  HO 
DR:  Range  increment  in  metres 

K0:  Final  value  of  height  in  ray  incr«nent  in  metres 
The  ray  trace  equations  that  are  used  in  calculating  DR  for  the  various 
values  of  H  (I)  or  the  final  value  of  height  H0  at  the  observing  range  are  given 
along  with  illustrations  in  Figure  B-6.  The  overall  algorithm  that  increments 
AS  and  computes  the  ray  trajectories  is  best  described  by  the  program  listing 
itself. 

Reference 

1.  Thompson,  William  I.,  Ill,  "Atmospheric  Transmission  Handbook,"  Transportation 
Systems  Center  DOT-TSC-NASA-71-6,  February  1971. 
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Case  I I  Ray  penetrating  linear  region  (upgoing  or  downgoing). 


H(l) 


H(I+1) 


/  2  -6 

A1  =IaO  +  2  X  10  (RM(I)-RM(H-l)) 

liK  =  Ul  -  A0)*Ra(l) 


f 


Case  III  Ray  maximizing  (or  minimizing)  within  linear  region. 


Al  s  0 

2 

=  H(I+1)  -  AO  ‘HSO/a 
DR  s  •  RS(I) 


Case  Ill!  Known  OR  but  unknown  (upgoing  or  dovngoing). 


Ai  =  AO  +  Da/as(i) 

2  2 

Hf  =  K(I+1)  +  R3(I).(A1  -  AO  J/2 


Figure  B-6.  Formulas  and  Illustrations  for  the  Three  Cases  that  Need  to  be 
Considered  in  Calculating  a  Bay  Trajectory. 


PR06PAH  image 

COdMON/OUM/  ISHIPCGiZO) ,LABFL<16) 

OIMEMSTON  H(5),T(5>,RM<3) .RSCS) 

DATA  I  SH I P /n.  1  CO,  100, 200, 0, 1  00, 1  go,  2  02  .P,  100,  ICO.  204.7,100.^ 

1 1 00, 206. 6, 1  CO,  100, 208, 5, 1 00, 10  J,21 1,4, 100,1  no,  2U,  3, 100, 100,21  7, 
S2, 100, ICO, 2 21 ,1.1 00, 100, 22 6, 20 ,60, 120, 16 0,^0, 60, 120, 160, 

120, 60,  120,160, 20, 60, 120, 16i., 20, 60, 120, 160,20,60,120, 140, 

120, 6G;  120. 140, 20, 60. 120. 14 0,20 .60,120,140.2  0.60,120,14  0/ 

DATA  LAPEL/'H*,'E',^I',*G' .  ', 

CALL  IN1TT(48C) 

5  TYPE  'ENTER  HEIGHT  <M),  TEMP(C),  DESCENDING  TO  ZERO — 5  LEVELS  MAX' 
DO  10  1=1,5 
ACCEPT  ,H<1) ,T<1) 

IF  CH(I  )  ,LT,  .001)  GO  TO  20 
10  CONTINUE 

CALI  NEWPAGE 
GO  TO  5 
20  IMAy=I 

?5  ACCEPT  'OPTICAL  WAVELENGTH  (MICRONS)  »  ',W 
ACCEPT  'OBSERVING  HEIGHT  (M)  »  '.MO 
ACCEPT  'RANGE  TO  OBJECT  Ufl)  =  ',PXM 
R=1.E7*RKM 
00  30  l=1,lMAy 
P*1COO.-.1  2*H(1) 

T<»?7’5  .♦Td  ) 

30  Rm(I)»77.6*P/TX*.564*P  /  (TX  <>U**2)  +  .157AH(  1) 

AS=-2*H0/P 
OELA  =  C.5/fi^ 

IMAXl  =  IMAX-1 
DO  40  1  =  1  ,  IM  AX1 

40  RS{I)=1 .E46*<H(I)-m(I*1))/ <RM(I)-RM(14l)) 

DO  50  1  =  1,  IMAXi 
IF (HO .GT.H  (I> 1) )  60  TO  60 
50  CONTINUE 

6  0  RMOsRM< 1,1) +1 ,E  +  6*(H0-H(I*  1))/RS(I) 

10*1 

CALL  NEWPAGE 
CALL  MOVARS  (90,200) 

CALL  MOVARS  (90,200) 

CALL  DRWARS  (100.200) 

DO  300  1*1.5 
CALL  DRWREL (0,100) 

CALL  DRWREL  (-10.0) 

30r.  CALL  DRWREL  (10,0) 

CALL  MOVARS  (1  OC, 190) 
call  ORwARS  (100, 200) 

DO  -^lO  1  =  1  ,  5 
CALL  DKWRFL(6  0,0) 
call  DRWREL  (C  .-1L) 

31D  CALL  DRWREL  (0,10) 

Ihmax  =  2'' 

I F (H ( 1 )  .GT. 25 . )  IHMAX=5n 
IF  (H  (1)  ,6T.  50.  )  1MMAX=10C 
YFACT=5C0./  IHMAX 

00  312  1=1.6 

IH»(I-1  )*1HMAX/5  .  ^ 

CALL  MOVARS  (45, 19P^(I-1)'1u0) 

call  ANMOOE 
312  •RlTt(1C,314>  IH 
314  FCR*"AT(  1X  ,  1?) 

DO  316  1  =  1, Ifc 

C»LL  MOVARS  <1  n  ,6?0-2U*l  ) 

CALL  ANMOOE 

316  .rfRITE(  ir,31P)  LAREL(I) 

31'-  FOR*"  AT  (  lx  .  A  1) 

R''MIN  =  26L/  . 

RMf(iAX  =  2P5  . 

XFACT  =  3CC./  (RMHAX-HmhIN) 

1X»100*  (RM<  II-RMMinI^XFACT  +  .S 
1  Y»?on+H(1 )  *YFACT-»  ,5 


BIO 


CALL  MOVABS (IX *1Y) 

CALL  WOVREL(>5*S) 

CALL  AK(«>ODE 

TYPF'w' 

CALL  *»OVAPS  (IX  .lY) 

00  320  I»2,  IMAX 
IXaipO*(RB(  I)-RMKlN)*XFACT+.5 
IY«200*H(I) *YFACT+.5 
320  CALL  0RWA8S  <1X,IY> 

IXsIOO^Td  )  *12*,5 
IY»200*H(1) *YFACT*.5 
CALL  MOVABS (1X-S,1Y^5) 
call  ANMODE 
TYPE  'T' 

CALL  MOVABS CIX ,1Y) 

^>0  322  t  I'*''* 

IX*100+T(I) *12*,5 
IY«200*M(I>*YFACT+.5 
322  CALL  DRUABS (IX ,2Y) 

CALL  MOVABS<0«UC) 

CALL  ANMODE 

TYPE  '  HEI6HT(M)  TEMP(C)  M  UNITS' 

DO  330  1  =  1 ,  IMAX 

330  MRITE(ir,3AC)  H < I > ,T ( I ) , RM (l ) 

340  F0BMAT<5X,F4.1.6X«F5.1.6X«  F5.1) 

CALL  MOVABS (500,140) 
call  ANMODE 
,,,  i.RlTE(ig,342)  M 

342  FORMATC*  optical  WAVELENGTH  ■',F5,2,'  MICRONS') 
call  MOVABS (500.110) 
call  ANMODE 
WRITe<1C,344)  HO 

344  FpRMATC'  OBSERVING  HEIGHT  ',f5.1,'  METERS') 
call  MOVABS (500,80) 

CALL  ANMODE 
-*RITE<ig,34<^)  RXM 

346  FORMAT('  RANGE  TO  OBJECT  ',F5.1,'  AM') 

CALL  MOVABS (770,220) 

CALL  ANMODE 
type  'OBJECT' 

IXR*700 

1YR«26Q 

DO  35P  IS*1  .20 
IH«2*IS 

CALL  M0VABS(IXR*ISHIP(1 ,I S) ,I YR+IH) 

CALL  0RwAE'S(IXR^ISHIP(2.IS),IYR^Im) 

CALL  «CVABS(IXR*ISHXP(3,IS) .lYR^IH) 

CALL  ANMODE 
TVPE  'IMAGE' 

I  rR«46G 
I  H»0 

BEGIN  PAY  TRACE 
7'J  AOsAE 
RTOT=n, 

1  =  10 
HC  =  HO 
PMn*RHr! 

':o  1 F  ( At'i.LT.n. )  60  TO  110 

TFST=A0**2*?,E-6*(HM(i)-rm0) 

TF  (TFST.LT.n.)  GO  TO  90 

A1=S0RT(TEST) 

tRr  (AI-AO)  *RS  <1  ) 

IF  (RTOT+OR  .GT  .R  )  GO  TO  20  0 
AOsAI 
Hfl  =  H  (  I  ) 

RTOT»PTOT+DP 

BHOsijH  (I  ) 

IFd.ECi.l)  60  TO  ?00 
1  =  1-1 
GO  TO  8C 
9(1  DR»-A0*RS  ( I  ) 


IF(RTOT^OR.6T.R)  60  TO  2Q0 
Mn«MO-Ar**2*RS<I>/2 
R»O«RF"C-A0**2/2.E-6 
RTQTsRTOT+DR 

11U  TPST  *  AC**?+2.E-A*<R‘*(If1  )-R»*0> 
lF(TESTrLT  .0.)  60  TO  120 
Ala-SQRT(TEST) 

0R*( A1-A0) *RS  <I > 

IF(RTOT+DR  .6T.R)  GO  TO  200 
PTOT*RTOT4DR 
1*1  +  1 

IFd.EO.lPIAX.AND.lH.EO.O)  GO  TO  22S 
IFd.EQ.lFlAX)  GO  TO  ?20 
A0*A1 
HQ*H( I > 

RM0=RM (I) 

60  TO  110 
120  0«>*-AO*9S(I  ) 

IF(RTOT+OR  ,GT.R)  60  TO  200 
hO*hO-AO**2*PS< 11/2. 
Rf'0*R«0-A0**2/2.E-6 
RTOT«RTOT+OR 

a0*0. 

60  TO  8C 

200  A1«A0+(R-RT0T)/RS(I) 
ftM*( A1 **2- A0**2 )/2. 

HO*HO+RSf  I  )*ONl 
IF(H0*6T.2G  •)  GO  TO  220 
IS*hO+1 

CAt.L  KOVAaSdXR+ISHloddSltlYR+lH) 
CALL  0R«ABSdXR*ISHlP(3d  S)  ♦lYR  +  lH) 
CALL  t«rvAeSdXR  +  ISH2P(3.I  S).I  YR  +  IH> 
CALL  0RwABSdXR4lSHlP(A»l  S)  dYR  +  IH) 
220  lHaIH+1 
225  AS«AS+OFLA 

IF(AS.LT.3  .E-2)  60  TO  70 
CALL  TIMPUT  <lCMAf;) 

call  NEWPAG 

IFdCMAR.Efi  .26)  60  TO  25 
60  TO  5 
EAfO 
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APPENDIX  C 

CSMPOTBR  PP06EAM  BAYPPOP 


This  appendix  descrilies  a  BASIC  computer  program  which  determines  whether 
a  refracted  line  of  sight:  path  exists  between  a  source  and  receiver  located  at 
arbitrary  heights  above  the  sea  surface. 

It  assumes  the  existance  of  a  oonst2uit  gradient  of  index  of  refraction  within 
the  reuige  ^  >  dn/dh  >  -1/a  (a  *=  radius  of  earth  =  6378.388  km) .  The  user  may 
either  Input  the  gradient  of  his  choice  of  default  to  the  value  for  the  standard 
atmosphere. 

As  shown  In  Figure  C-1,  the  following  outputs  are  available: 

1)  Distance  to  source  horizon  (km) 

2}  Dip  angle  to  source  horizon  <mrad) 

3}  Distance  to  receiver  horizon  (km) 

4)  Dip  angle  to  receiver  horizon  (mrad) 

5)  Beam  pointing  and  receiving  angles  for  the  refracted  line  of  sight 
path  (mrad) 

6)  Baytrace  from  source  for  arbitrary  beam  pointing  angle  and  divergence 

7)  Raytrace  printout 


Cl 


HORIZON 


Figure  C-1.  ScUi^le  Outputs  and  Raytrace  for  the  RAYPROP  program. 


4  RBY  •**RATPPOP** 

5  INIT 

6  PAGE 

7  GO  TC  10a  -  - - - - -  - 

S  LET  G=-2.3S“e 

f>  GO  TO  2^0  - - -  ...  - 

1?  GO  TO  363 

IP,  50  TO 

2?  LET  I3«l 

21  17  ■P>2  THEK  553 

2?  GO  TO  413 

'  •100  HEM  ■  ...... 

102  LFT  P=3 

110  SET  KEY  - 

1?3  PPrn  "ITPUT  SOURCE  AND  RECEIVER  HEIGHTS  IS  METERS  " ; 
130  INPUT  H0,H1 

140  PRINT  *’j  INPUT  PATH  LENGTH  IB  KILOMETERS  • 

150  INPUT  §  ”  . . 

163  PRINT  ‘j  INPUT  INDEX  GRADIENT  PER  METER.  OR  PRESS" 

161  ttRINT  "USER  KEY  2  FOR  STANDARD  AIMOSPHERE 
173  prv  USE  NEGATIYt;  GRADIENTS  FOR  INDICES 
183  REM  WHICH  DECRE.ASS  WITH  HEIGHT. 

193  INPUT  G 

233  LET  A-6-576.36S 

213  LET  K«1/(1’^A*G*1330) 

215  17  K>3  ’^EEN  223 

216  PRINT  "j  " 

217  GO  TC  320  . - . — 

223  RIW  COMPUTE  H0RI20NS 

233  LET  Z3«SCR(A*K3*:/5-33) 
py'*  let  Zl»SCR(A*Hl*K/533) 

241  LET  30— Ea/Z3-530*70/(K*A) 

242  LET  1,1— H1/Z1-503*Z1/(':*A) 

253  WRINT  ",T  J  J  SOURCE  HCRISCK  *  ";Z0;"  KM' 

251  PRINT  , .J”SOURCE  HORIZON  ANGLE  »  ’;30;"  MRADJ  J  J  " 

260  PRINT  J”RECEIVER  HORIZOIL-*..  ^214  KM  ...  " _ .7.  _ 

?«1  PRINT  "J"RECEIVSR  HORIZON  AHGLE  »  "jBl;"  MRADJ  J  J  " 

270  REM  TEST  TOR  REFRACTED  LINE  OF  SIGHT  -  "  . 

280  IF  Z3+Z1>S  THEN  313 

293  PRINT  ’NO  RETRACTED  LINE  OF  SIGHT  PATH  EXISTSJ  J  J  " 

.330  GO  TO  A4? 

313  PRINT  "refracted  LINE  0?  SIGHT  PATH  EXISTSJ  " 

323  LET  I3={H1-E’>)/S-7.30*S/{K*A) 

321  LET  B1«(H0-H1 )/S-500*S/(K*A) 

333  PRINT  "beam  POINTING  ANGLE  =  MR.ADJ  " 

3.31  PRINT  "RECEIVING  ANGLE  =  ";B1:  MRADJ  J  J“" 

340  PRINT  PRESS  USER  KEY  3  FOR  A  RAYTRACE  0R~KEY  4  TO  EHDJ  J 
:  .  350  WATT  -  ..  -  ■  .  " 
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363  . MYTR^CE  ROUTINE 

376  PRINT  ’’input  BB4M  INCLINATION  ANT)  DI\ER&ENC?  IN  MRAD  ; 

380  HEM  ANGLES  ABOVE  HORIZONTAL  ARE  PCS  IT  IV I , BELOtf  ABE  NEGATIVE 

390  INPUT  Ba.D 

466  Jim  DRAW  SEA  SURFACE 

.  .465  LIT  13-32  - - - - - 

416  PAGE  616: 

426  REM 

430  WINDOW  0,S+2.0,12 

446  MOVE  6I6«1,0  . . .  ■  . . — ■  . . 

453  FOR  L-3  TO  S  STEP  S/200  . 

•  466  LET  X-L*l  >  ■  ... 

470  LET  T«500*L*(S-LVA 
486  DRAW  6I6:X,Y 
490  NEXT  L 

600  REM  DRAW  SOURCE  AND  RECEIVER 

.513  MOVE  410;  1.0 

523  DRAW  010:1, H3  -  - - - 

533  ROVE  4I0:l+S,0 

546  DRAW  0I6:1+S,HI 

553  REM  PLOT  RATS 

566  FOR  B«36“D/2  TO  B3+D/2  STBP  D/2 

670  MOVE  0ie:l,R0 

560  FOR  L«0  TO  S  STEP  S/200  .  -  - 

593  LET  H=H0+3*L+500*L'‘2/(’"*.A) 

'  591  IE  H<0  THEN  655  - 

600  LET  x=L+l 

•  610  LET  Y-H+533*L*(S-L)/A  . . 

I  611  IF  T>10  TREK  655 

»  626  DRAW  fI0!X.F - - - . . . . . . . 

630  NEXT  L 

636  IF  D»0  THEN  636  •  -  - - - 

646  /lEXT  B 
653  HOME 
655  LET  P=f+1 

660  PRINT  -PRINTOUT-REY  5  -  ^  - 

99B  PRINT  'g  G  G  G  G 
999  END 
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